Frequency splitting of intervalley phonons in graphene 
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We study the thermal distribution of intervalley phonons in a graphene sheet. These phonons 
have two components with the same frequency. The degeneracy of the two modes is preserved by 
weak electron-phonon coupling. A sufficiently strong electron-phonon coupling, however, can result 
in a splitting into an optical and an acoustic phonon branch, which creates a fluctuating gap in 
the electronic spectrum. We describe these effects by treating the phonon distribution within a 
saddle-point approximation. Fluctuations around the saddle point indicate a Berezinskii-Kosterlitz- 
Thouless transition of the acoustic branch. This transition might be observable in the polarization 
of Raman scattered light. 
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Graphene, a two-dimensional sheet of carbon 
atoms forming a honeycomb lattice, has highly un- 
usual electronic properties [ll-13| . This is due to the 
fact that there are two bands that touch each other 
at two Dirac nodes. Moreover, the low-energy 
quasiparticles of undoped graphene experience a 
linear dispersion around the Dirac nodes. Trans- 
port properties, characterized by the longitudinal 
conductivity at the Dirac nodes, are quite robust 
and do not vary much from sample to sample. Ex- 
actly at the Dirac point a minimal DC conductivity 
has been observed in a number of experiments llj- 
1^ . Its value is not much affected by disorder [s!, ^ . 
A typical source of disorder in graphene are frozen 
lattice deformations (ripples) |3, §]■ In a more re- 
alistic description, these deformations may not be 
frozen but fluctuate thermally due to the softness 
of the two-dimensional material. 

In this paper we investigate effects of the 
electron-phonon interaction in graphene. Con- 
sidering only in-plane displacements of graphene 
atoms, we have four different phonon modes. For 
low-energy electronic states electron-phonon inter- 
action is efficient if the phonon wavevector is close 
to r, K or i^T'points. As it was explained by Basko 
and Aleiner [9|, the only modes which are effec- 
tively coupled to electrons are the pseudovector 
optical phonons corresponding to the irreducible 
representation E2 of the group C^y from the F 
point, and the scalar phonons corresponding to the 
irreducible representations Ai and Bi of the same 
group from the points K and K' . 

Experimental information on the interaction is 
obtained by Raman spectroscopy [lol - [l^ an d by 
angle-resolved photoemission spectroscopy [l5|. 
The above mentioned theoretical analysis of the 
role of different phonons is confirmed by the fact 
that in the Raman spectrum of graphene only two 
two-phonon peaks are seen: the D* peak and the 
G* peak, corresponding to scalar and pseudovector 



phonons, respectively. In our paper we will con- 
sider interaction of electrons only with the scalar 
phonons. 

The physics of the electron-phonon system is de- 
fined by the Hamiltonian with optical (Einstein) 
phonons at frequency ujq [16|: 



r')]4cr 



(1) 

Here cj (cr) are the electron creation (annihila- 
tion) operators and bl ^, (6r,r') are the phonon 
creation (annihilation) operators. This Hamilto- 
nian describes an effective attraction between the 
fermions and a renormalization of the electronic 
hopping amplitude /ir,r', reducing the hopping rate 
substantially (polaron effect). The attractive inter- 
action may lead to the formation of Cooper pairs. 
In 2D this effect is not relevant due to strong fluc- 
tuations, preventing the system to become super- 
conducting. 

The conventional approach to determine the 
properties of the phonons and the electrons is 
based on a self-consistent evaluation of the self 
energy (Migdal approximation) [13, [13 ■ The lat- 
ter provides an effective (or renormalized) en- 
ergy and its imaginary part an effective scattering 
rate. Such a static approximation might be insuf- 
ficient in a two-dimensional system, since it does 
not take into account thermal fluctuations. This 
was already discussed in an experimental study 
of graphene (Tsj . To avoid this problem, we in- 
clude thermal fluctuations in our approach. To this 
end, we replace the phonon operators &r,r', bl ^, by 
their quantum average: 6r,r' ~ {br.r') = Vr,r' and 

r' ~ i^l r') = ''^r r' • ^^^^ approximation we 
can keep thermal fluctuations but ignore quantum 
fiuctuations of the phonons. The electrons, on the 
other hand, are studied in full quantum dynam- 



2 



ics. This reduces the grand-canonical ensemble at 
inverse temperature /3, defined by the generating 
function Tre~^^ to a functional integral with re- 
spect to thermal fluctuations of the lattice distor- 
tations Ur,r' and a trace with respect to the elec- 
trons. It should be noticed that only the real part 
Re{u) couples to the electrons. After performing 
the trace over the electrons we get 



with 



det(l + e-^'')e-^^°V[u] = Z. (2) 



with the dispersion for the rescaled phonon field 

C _ "^0 \ " 2 
Do - ^ 2^ U^y . 

r.r' 

It is convenient to introduce a sublattice rep- 
resentation for the tight-binding Hamiltonian h 
and to expand it around the two valleys K and 
K': The graphene unit cell contains two atoms, 
each of them has one 7r-orbital. This gives a two- 
component wavefunction. Moreover, the band- 
structure has two nodes (or valleys). Expansion 
around the valleys leads to a wavefunction that 
is represented by a 4-component column vector. 
We shall work in a direct product vector space 
of the valley space and the sublattice space. For 
the phonons we consider here only the interval- 
ley contribution, which play a major role in the 
electron-phonon interaction This gives for the 
low-energy Hamiltonian 



u ■ nSo 



(3) 



The parameter t is the bandwidth. Here we have 
used a coordinate system that refers to one sub- 
lattice and one valley in the notation of Ref. fio'l : 
the indices of the Pauli matrices Ej (j = 0,...,3) 
are acting on the A-B sublattices and the indices 
of the Pauli matrices {lij {j ~ 0, ...,3) are acting 
on the two valleys K and K'). Moreover, we have 

$ — Tj ■ where 9 is a lattice difference operator. 
For graphene in a homogeneous magnetic field this 
electron-phonon interaction was considered in Ref. 
[20j . Then the intervalley phonon field has two 
component: 

^/ ) ^ "l^rlli + U2,rll2 , 

where uir — {uy. + u'^)/2 and uir — ~i{ur — u'j.)/2 

The integral Z in Eq. ([2]) serves as a generating 
function that allows us to get, for instance, the 
static electronic Green's by differentiation of In Z 
as 

Gr,r' = | j{l + e-P%^^,e-0SV[u] (4) 



5 = 5o-/3-Mogdet(l + e-'3''), u^im^u^). (5) 



It is useful to notice that e 



-f^^o detfl -h 



g-Ph^ is a non- negative function. Therefore, 
g-0S 1^ is a probability density for the phonon 
field, and the static one-particle Green's function 
then can also be written as an average {■.■)ph with 
respect to the distribution e~^^ /Z [2l| : 



Gr,r'=((l + e-^")-J,U 



(6) 



The distribution e~^^ jZ is invariant under a 
unitary transformation. For the special transfor- 
mation 



(7) 



with < < 27r the electron-phonon coupling 
term in Eq. (jSj satisfies the relation 



Uu ■ mj'^ = Om • H 



(8) 



where O is the orthogonal transformation (i.e. ro- 
tation by angle (p) in the space of the matrices H: 

O = cos(y5HoSo + isinipIl2T,Q — exp(i(y5H2)So ■ 

Since the gradient part of the Hamiltonian h is 
invariant under a global rotation of the phonon 
field u, the distribution e~^^ jZ is also invariant. 
As a consequence, the phonon field can produce 
massless fluctuations. 

Second order perturbation theory is the stan- 
dard approach for evaluating the change of the 
phonon frequency by the electron-phonon interac- 
tion [7, 22, 23]. In our model, defined by the action 
([S]) , second order in a gives us for the renormalized 
frequency of the intervalley mode 



M=0 



(wo - a2)8v 



with 



a 



In 



coshifiEF) + cosh(;3tA) 
1 + cosh{PEF) 



(9) 



(10) 



where A = 2^/^T is the momentum cutoff. For fit ^ 
1 , which is satisfied even at room temperature due 
to i ~ 2.7 eV, the renormalized phonon frequency 
has the asymptotic behavior 



Thus the intervalley phonons are softened by the 
electron-phonon coupling. Moreover, the phonons 
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hardens as we go away from the Dirac point 
Ep = 0, in agreement with recent experiments 
[3 113, Hi • However, this frequency becomes neg- 
ative if Wo < 02, indicating an instability of the 
electron-phonon system. In particular, the vanish- 
ing frequency reveals a phonon softening, where 
the optical phonon mode becomes acoustic. In the 
following we shall discuss that this instability is as- 
sociated with a splitting of the degenerate phonon 
modes, where instead of a single phonon frequency 
two different frequencies appear. 

The instability cannot be treated within second- 
order perturbation theory but requires a self- 
consistent approach. Here it is natural to perform 
the integration in Eqs. ^ and (|H) in a saddle- 
point (SP) approximation. This leads to the SP 
equation SuS = which determines an average 
lattice distortion \u\ — u. Assuming a uniform 
solution u the SP equation dS/dui — reads 



2 — u 

UJQ 



sinh(/3efc) 



Efc cosh(^£;F) 



-cosh(/3efe) (2^)2 

(11) 



with the dispersion e/c — y/WW + ^P of the elec- 
trons. A nonzero u opens an electronic gap A = 
2u. This equation fixes only u, not the direction of 
the vector uq = (■ucos(/3, usimp). For a non-trivial 
solution u we we can perform the integration 
in Eq. ((TT]) and obtain the equation 



cosh(/?£;f ) + cosh(/3VA2^2 ^ ^2) ^ ^^^^^2^^/^2 
cosh^PEp) + cos'h{f3u) 

(12) 

Then the critical point Pc of the instability is de- 
termined from this equation for u ~ 0. 
For /3t > 1 and \Ef\ < m we get from Eq. 



2a2 / 7rt2^2 

U ~ 1 



(13) 



and for f3t > 1, \Ef\ > u 




(14) 



While an increase of the electron-phonon coupling 
a increases u for \Ef\ < S, it decreases u for 
\Ep\ > u. This result indicates that the insta- 
bility of the u — solution is supported by the 
electron-phonon interaction near the Dirac point, 
whereas the instability is suppressed further away 
from the Dirac point. 

For small u Eq. (fT2|) can be expanded. This 
allows us to study the properties of the phonon 
renormalization for nonzero u, i.e. in the regime 
where u — and the perturbation theory are un- 
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FIG. 1: Splitting of the plionon modes (from Eqs. (|17l 
((TH|). The red (blue) curve represents lji (cj2) in units 
of iJo, indicating that for the dimensionless electron- 
phonon coupling q/ V tiJo > Oc ~ 0.94 there exist two 
different frequencies. 



stable. The SP equation then reads 
o dS 



(w — a2 + a^u )uj — . (15) 



Besides the solution u = this equation has the 
nonzero solutions 



±y/a2 - Wo/ (a2 - Wo > 0) . 



with 



duiduj 



4(a2 - wo)- 



(16) 



instead of Since the SP solution uo is on a 
circle with radius u, the eigenvalues of the matrix 
on the right-hand side are wi — 4(a2 — wo) and 
W2 = 0. For a stable solution, on the other hand, 
u = the frequency shift agrees with the result of 
the second-order perturbation theory. Combining 
(O (for 02 < Wo) and (for 02 > Wo), we obtain 
the renormalized phonon frequencies as 

w |a2 - wo| [1 + 6(02 - Wo)] , (17) 



W2 ~ (wo - a2)6(wo - 02) . (18) 

Thus, the mode with W2 vanishes for the solution 
w 7^ 0, as plotted in Fig. [TJ 

After having solved the SP equation for a uni- 
form phonon field u, we study the fluctuations 
around the uniform solution. Fluctuations may 
play an important role and should modify the 



4 



uniform solution due to the two-dimensionality of 
graphene. In particular, the fluctuations are cru- 
cial for the vanishing phonon frequency uj2 in Eq. 

m- 

With the special SP solution uq also any rotated 
phonon mode Ouq is a SP solution. Therefore, we 
can choose the special solution uq = (u, 0) and, 
using the identity ([8]), obtain the action at the SP 
as 



5.5o-ilo, 



det(l 



im3$+uU'aiU'' 



which is degenerated with respect to a global uni- 
tary transformation U . In order to study fluctua- 
tions around the SP solution we introduce a spa- 
tially fluctuating unitary matrix Ur- This can be 
cast into a nonlinear matrix field 



$r = C/rniC/t 







(19) 



such that the action becomes 



5: 



5o-^log 



det(l 



-0(-itn3^+G*)^ 



) ■ (20) 



This action can be expanded, either for u/t < 1 or 
for u/t > 1. 

Now we assume that u > t. Formally, t can be 
chosen independently of the real bandwidth of the 
system. But this choice means that the physics is 
restricted to quasiparticles up to energy t. Then 
the expansion of the action S of Eq. (|20l) in powers 
of 1/u provides a C/(l) nonlinear sigma model: 



S^S'„+ F{P, u)Tr 
with the prefactor 



(21) 



u 1 , , 
- + - n(l 



4'/?-- "21.^^ 

(22) 

Using the result of u for /3m > /St ^ 1, the (3 
dependence drops out of the prefactor: 



t^UJQ 



1 



(23) 



This result can be considered as a renormalization 
effect for the temperature, where the renormalized 
dimensionless temperature reads 



u 



2^2 



Wo 



i2 



1 - 



2, ,2 



4a4 



T 



(24) 



Here we have assumed |i?_F| < u and have used 
the expression of u in Eq. ([T^. As the electron- 
phonon coupling a increases, the renormalized 
temperature r increases as well. This reflects the 



fact that the electronic fluctuations enhance the 
phonon fluctuations. The case \Ep\ > u is rather 
unrealistic here, since we also have assumed u > t. 

The fluctuating term in Eq. ([21]) represents an 
XY model for the angular fluctuations ipr- Thus, 
the phonon fluctuations undergo a Berezinskii- 
Kosterlitz-Thouless (BKT) transition if Tc is of 
order one [25| . Thus, for temperatures below 
the BKT transition point Tc the fluctuations are 
strongly correlated whereas above this tempera- 
ture the correlations of the fluctuations decay ex- 
ponentially due to the proliferation of vortex pairs. 

We conclude from our calculation that second- 
order perturbation theory with respect to electron- 
phonon interaction reveals a phonon softening. 
This result is in agreement with other calculations 
0, m, H^l • -fo'" sufficiently large electron-phonon 
coupling the perturbative approach breaks down 
and a self-consistent calculation is necessary. We 
have used an SP approximation for the phonon 
fluctuations and found an instability of the per- 
turbative approach due to a splitting of the opti- 
cal phonon mode into an optical and an acoustic 
branch. The latter is related to a massless mode 
caused by the rotational symmetry of the two sys- 
tem. It is characterized by long-ranged correla- 
tions of the phonon fluctuations, in contrast to 
the short-range fluctuations of the optical phonons. 
However, at sufficiently high temperatures these 
long-range correlations can undergo a BKT transi- 
tion by the creation of vortex pairs, resulting again 
in short-range correlated fluctuations. The aver- 
age lattice distortion vanishes for all regimes as a 
consequence of the Mermin- Wagner argument. A 
possible way to observe the BKT transition exper- 
imentally is to measure the polarization in Raman 
scattering, since the photon polarization couples 
to the direction of the lattice distortions. 

The instability of the intervalley phonons is 
very different from the instability of the Holstein 
(out-of-plane) phonons studied recentl y, w here the 
phonons undergo an Ising transition [26|. In the 
latter case only short-range correlated fluctuations 
appear (except for the critical point) but the aver- 
age distortion has a non-zero value. 

To understand the effect of the thermal phonon 
fluctuations on the electronic transport properties, 
we can return to Eq. ^ . Then produces a ran- 
dom gap in the electronic spectrum which has to 
be averaged with respect to the phonon distribu- 
tion e^^^ /Z. Thus, the effect of thermal phonons 
is similar to the effect of frozen correlated disor- 
der, depending on the temperature of the sample 
though. Since an uncorrelated random gap with 
vanishing mean does not affect the minimal con- 
ductivity [27l - [29| . we expect a similar result for 
thermally fluctuating intervalley phonons, at least 
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near the Dirac point. 

In conclusion, our calculation reveals that the 
electron-phonon interaction in graphene leads to 
a substantial renormalization of the intervalley 
phonons. For sufhciently strong electron-phonon 
coupling this causes an instability, where the 
phonon frequency vanishes and a new pair of 
phonons appears, consisting of an optical and an 



acoustic branch. 
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